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Abstract. For integers n, q = 1,2,3, . . ., let Pol,,,, denote the C-linear space of 
polynomials in z and z, of degree < n - 1 in z and of degree < q - 1 in z. We 
supply Pol„ , ; with the inner product structure of 

^ L 2 (c,e~"' |z|2 cL4), where dA(z) = z = x + iy; 

the resulting Hilbert space is denoted by Pol m ,„ i9 . Here, it is assumed that m 
is a positive real. We let K m ,n,< ( denote the reproducing kernel of Pol ra ,„^, and 
study the associated determinantal process, in the limit as m, n —> +oo while 
ly-j n = m + O(l); the number q, the degree of polyanalyticity is kept fixed. We 

call these processes polyanalytic Ginibre ensembles, because they generalize 
the Ginibre ensemble - the eigenvalue process of random (normal) matrices 
with Gaussian weight. A possible interpretation is that we permit a few higher 
Landau levels. We consider local blow-ups of the polyanalytic Ginibre ensem- 
bles around points in the spectral droplet, which is here the closed unit disk 
D := jz e C : |z| < 1). We obtain asymptotics for the blow-up process, using a 
blow-up to characteristic distance m~ 1/2 ; the typical distance is the same both for 

S interior and for boundary points of D. This amounts to obtaining the asymp- 
totical behavior of the generating kernel K,„^ q . Following |2|, the asymptotics 
of the K m/n< o is rather conveniently expressed in terms of the Berezin measure 
I (and density) 

AB% hq (w) := B«^(B7)dA(H»), B«» = 'W^ 1 e -»tf 

P^m,n,q\Z, Z) 

Q\ For interior points |z| < 1, we obtain that dB^ ^w) — > d<5 2 in the weak-star sense, 

where <5 Z denotes the unit point mass at z. Moreover, if we blow up to the scale 
of m~ 1/2 around z, we get convergence to a measure which is Gaussian for q = 1, 
but exhibits more complicated Fresnel zone behavior for q > 1. In contrast, for 

t-H exterior points |z| > 1, we have instead that dB^ ? (io) — > dw(w, z, D 1 '), where 

dw(zv, z, D 1 ') is the harmonic measure at z with respect to the exterior disk 
V := \w £ C : \iv\ > 1). For boundary points, |z| = 1, the Berezin measure dB^ > H1/ 
converges to the unit point mass at z, like for interior points, but the blow-up to 
the scale m~ 1/2 exhibits quite different behavior at boundary points compared 
with interior points. The Fresnel-type pattern appears also for boundary points 
when q > 1, but then it is not rotationally symmetric. 
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1. Introduction 

Notation. We will write use standard notation, such as dX and int(X) for the 
boundary and the interior of a subset X of the complex plane C. The complex 
conjugate of a complex number z is usually written as z. We write 1R for the real 
line, D := \z 6 C : |z| < 1} for the open unit disk, and D e := {z e C : |z| > 1} for the 
open exterior (punctured) disk. The characteristic function of a set E is written 1 E . 
We write 

dA(z) = 7i _1 dxdi/, where z = x + iy e C, 
for the normalized area measure in C, and use the standard Wirtinger derivatives 

d z := ~(d x - id y ), d z := ~(d x + \d y ), where z = x + iy. 
We also write A for the (quarter) Laplacian 

A 2 :=d z 3 z = \(d 2 x + d)). 

Determinantal projection processes. Given a locally compact topological space 
X with a Radon measure \i, a determinantal projection process (in the sequel 
just determinantal process) is a random configuration of n points defined by the 
following probability measure on X": 

(1.1) dP(zi,...,z n ) = ^det[K n (z tr z j )]l j=1 dfi(z 1 )---d^(z n ). 

Here, K n is the integral kernel of a projection operator to an n-dimensional sub- 
space of L 2 (X, f.i). It is customary to identify all the permutations of the points and 
think the process as a random measure X!/=i ^z, on X. 

A general definition of a determinantal process was introduced by Macchi 
|1T7| , who wanted to model fermions in quantum mechanics. Indeed, for any 
determinantal process, the probability density vanishes whenever any two points 
in the n-tuple (z\, ...,z„) coincide (fermions are forbidden to be in the same state). 
We interpret this as saying that the points in the n-tuple repel each other. Point 
processes of this kind appear in several contexts, e.g., in random matrix theory 
and combinatorics (for general surveys, s |15], |8|; we should also mention the 
books flg|, 0, GUI, El). 

Eigenvalues of random normal matrix ensembles. Our main motivating exam- 
ple comes from the theory of random normal matrices. This topic has in re- 
cent years been subject to rather active investigation by physicists as well as by 
mathematicians. For an introduction, see, e.g., fl24ll . So, we shall use X = C and 
d/./(z) = e~ m ^ z 'dA(z), for a positive weight function Q satisfying some mild regu- 
larity and growth conditions; m is a positive real parameter, and cL4(z) = n~ l dxdy 
is the normalized area measure. Let us write L 2 (C,e _mQ ) := L 2 (X,,u) in this situ- 
ation. The determinantal projection process is associated with an n-dimensional 
subspace of L 2 (C, e~ m< ^), and we will use the space Pol„ of all polynomials in z 
of degree < n — 1; we write Pol„ M , to indicate that we have supplied Pol„ with 
the Hilbert space structure of L 2 (C, e - '"^). The density of the process is then given 
by the reproducing kernel K„ hn of the space Pol„ M ,. So, we are talking about the 
probability measure 

(1.2) dP( Zl , . . .,z n ) = i det[JC ffl ,„(z„z ; )];.; y=1 e- m| Q( 2 ^ + ■ + Q( 2 »)»dA(z 1 ) • ■ ■ cL4(z„). 
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In terms of the correlation kernel 

(1.3) m min (z,w) := K m , n (z,zv)e- 1 2 m{Q(z)+Q(w)] , 

which is the an integral kernel of an orthogonal projection L 2 (C), the expression 
( |1.2) simplifies to 

(1.4) dP( Zl , . . . ,z„) = - det [yK m , n (z i/ Zj)] n i . =i dA(z 1 ) ■ ■ ■ dA(z„). 

The process described by ^1.2) and | |1.4[ represents the eigenvalues of a random 
normal matrix picked from the distribution 

(1.5) ^e- mtr[ Q( M ) 1 dvol nm()l) (M) / 

where dvol nm („)(M) is the natural Riemannian volume form on the n X n nor- 
mal matrices inherited from the metric of C" ; Z ni „ is the normalization constant 
needed to make the total mass equal to 1. We are interested in the limiting be- 
haviour of the process as m,n — > +oo while n = mx + O(l) for some positive real 
number t. Without loss of generality, we will consider only t = 1. 

Local blow-up processes. Let N+ and N+ t o be the set of points defined by 
N+ := {iveC: AQ(w) > o} ; N + , Q :={weC: AQ(iv) > o}. 

In the arXiv preprint Il2l , which will appear later in the expanded form 1 13 |, the 

function Q was defined as a certain envelope of Q, namely the largest subharmonic 
function in C which is < Q everywhere and has the growth bound 

Q(z) = log|z| 2 +0(l), aslzl^+cx,. 

It is known that AQ = l^AQ for some compact set S (see, e.g., [13J). We assume 
that S is the minimal compact with this property, and call S a spectral droplet. 
We then have S c A/V,o- The point process {1. 2) has the following property: as 
m,n — » +oo while n = m + O(l), the points will tend to accumulate on the set S 
with density AQ there. Moreover, the set S n N+ is rather regular for real-analytic 
Q, as the Sakai theory applies (cf. [14J). Typically we then expect a real-analytic 
boundary, with the exception of cusps and contact (or kissing) points. Let us refer 
to the set int(»S n W+) as the bulk. The results of |2], |3] show that for bulk points z, 
the local blow-up process at z, with coordinates . . . £n), 

^:=m 1/2 [AQ(z)] 1/2 (z i -z), 

where (z\, . . . ,z n ) are from the process \1.2\ , converges weakly to the translation 
invariant Ginibre(oo) process, as m, n — > +oo while n - m + o(l). The associated 
generating kernel is the reproducing kernel (£, rf) i-» of the Bargmann-Fock 
space. This has the flavor of a universality result. The corresponding statement in 
the Gaussian unitary ensemble (GUE) case is the universality of the sine kernel for 
bulk. We observe here that the sine kernel is the reproducing kernel for the Paley- 
Wiener space (a subspace of L 2 (R) consisting of entire functions). As for the two 
boundary points in the GUE model, the Tracy-Widom distribution appears, which 
is generated by the Airy kernel. The Airy kernel is reproducing for another Hilbert 
space of entire functions. This suggests that for real-analytic Q and z 6 dS D N+, 
there should exist a local blow-up 

£;:=m e (z,--z), 
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where 8 = 8(z) is a suitable positive real, such that as m,n — > +00 while n = 
m + o(l), the process (£1, . . . , £„) would converge to a determinantal process whose 
generating kernel is the reproducing kernel of a Hilbert space of entire functions. 
We verify this in the context of the Ginibre ensemble (i.e., with Q(z) = |z| 2 ), and 
identify the associated Hilbert space with a closed subspace of the Bargmann-Fock 
space characterized by slow growth in a half-plane. In that case, 6 = \ as in the 
case of interior points. 

The Berezin measure and the Berezin density. In |2J, [3J, Ameur, Hedenmalm, 
and Makarov study the Berezin measure dB^,, and Berezin density E^f n : 

jd(z) ( 7 ,,\ ._ R <z> ( 7n \AA( 7n \ R<z> ( 7n \ ._ l-KnV'C 2 ' W )\ 2 -mO(m) 

^m,n\Z, Z) 

which arise in the study of the Berezin transform. The Berezin measure is a 
probability measure, which makes it more stable than the reproducing kernel 
K mi „ of Pol m/ „ itself as we let m,n — » +00. In terms of the point process \l-2) , Ej£„ 
measures the amount of repulsion from z caused by placing one of the points at 
z. For bulk points z, we have the convergence dB„ „ — » dS 2 in the weak-star sense 
of measures, as m,n —> +00 while n - m + o(l). Here, db z is the Dirac point mass 
at z. In fact, there is a better result: the blow-up Berezin density 

B< z > (£) = I E< z> (z + g \ 

mAQ(z) \ y?wAQ(z) ' 

converges to the standard Gaussian e~'^' 2 . This corresponds to the convergence of 
the local blow-up of the point process to the Ginibre(oo) process (cf . [2 1, [3 ]). On the 
other hand, for points z outside the spectral droplet, i.e., for zeC\5, the Berezin 
measure dB^„ converges in the weak-star sense of measures to harmonic measure 
at z with respect to the exterior domain C \ S as m,n — > +00 while n - m + o(l) 
(see |2] for Q - |z| 2 , and |3J for the general result). 

The local blow-up of the point process and the Berezin density. It is convenient 
to think of the point process | |1.2[ in terms of the A:-point intensities 

det [>K m ,„(z,,z ; )])\ =i = det [^(^zyJe-HQfeHQ^)]]*^. 

We notice quickly that the intensities are unchanged if the kernel changed to 

Ki,„(z r w) := x{z)x{w)K mi „{z,w), 

provided that x is measurable with \x(z)\ = 1 (we can call this a "guage transfor- 
mation"). This can help in the asymptotical analysis of local blow-ups. For k - 2, 
we get the 2-point intensity 



[K m ,n(Zl,Zl)Km,n( z 2/ z 2) ~ l^m,n( z l/ z 2)| 2 }e 



-mQ(z 1 )-mQ(z 2 ) 

I K(zi.zi) ) 



K(zi,zi) /' 

where we recognize the Berezin density as a correction to the product of the two 
1 -point densities, where the 1 -point intensity is 

K nhn (z 1/Z ,)e- m ^\ 
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So, as far as the 2-point intensity goes, we just need the 1-point intensity and 
the Berezin density. Since the 1-point intensity is just the restriction to Z\ = 22 
of the Berezin density, the Berezin density is all we need to describe the 2-point 
intensity. We will be a little lazy and just work with the Berezin density in the 
context of local blow-ups, although the asymptotics of the fc-point intensity would 
strictly speaking require a little more work. So, although we state our many of our 
assertions regarding local blow-ups in terms of the Berezin density, we maintain 
that they generalize to statements about the point processes (cf. O). 

The Ginibre ensemble and its polyanalytic generalization. The case Q(z) = |z| 2 
of | |1,5) (or \1.2\ is known as the Ginibre ensemble. The (probability generating or 
reproducing) kernel is now particularly simple: 

K m ,„(z,ro) = m^_ t — ^ — • 
7=0 h 

Here, S = D, the closed unit disk. We will consider a family of generalizations of 
the Ginibre ensemble, the polyanalytic Ginibre ensembles, which are defined by the 
reproducing kernels K mAi(j of the subspaces 

Pol m< „, ? := span |z 7 z* : < / < tt — 1, <k < q -1^ 

supplied with the Hilbert space structure of L 2 (C,e~ m ' z ' 2 ). The parameter value 
q = 1 corresponds to the standard Ginibre process. In general, we now project to 
the n^-dimensional subspace of the polyanalytic polynomials, where the degree in z 
is < q - 1, and the degree in z is < n — 1. Note that the dimension of the subspace 
is now nq and not n. We can think of the case q > 1 as permitting more Landau 
levels than the lowest one (there are two similar models, see 1 5 1, 1 19 1, 1 25 ]; see also 
|1J). We will keep the polyanalyticity degree q fixed in the process, while we let 
both m, n tend to infinity. We remark here that we have noticed that Tomoyuki 
Shirai is interested in a similar model l22l . 

As for the point process, the points generally repel each other, but for q > 1, 
they also tend to avoid certain geometric configurations, such as circles and lines. 
We have run a simulation in Figure [TTT| 

Results. Macroscopically, we find that the behavior of the polyanalytic Ginibre 
ensemble is similar to that of the Ginibre ensemble. If we measure this in terms of 
the Berezin measure, we have that as m, n — > +00 while n = m + O(l), 

dBm,«,<? -> d6 z for z 6 D, dB%], hq da> z for z e D c , 

where co z is the harmonic measure of z with respect to the exterior disk D e . Here, 
K m ,n,q is the reproducing kernel for Pol„ wl/ ,j, and we use the notation 



*^m,n,q\Z, Z) 



for the corresponding Berezin measure and Berezin density. Interestingly, the 
microscopic behavior of the Berezin measure dB^ n ^ is quite different for q > 1 
compared with the Ginibre case q — 1. In terms of the blow-up Berezin density 
(AQ(z) = A|z| 2 = 1 here) 

BL z , > M (5) = m- 1 BLl ? (z + m- 1/2 5), 
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Figure 1.1. The polyanalytic Ginibre process with the kernel 
K m , n ,a with m = n = 61 and q = 3. The simulation is based on 
the algorithm by Hough, Krishnapur, Peres, and Virag |15|. 



we have the following asymptotics as m, n — > +oo while n — m + O(l): 

bL1,(«) — » ^L^fl^V 1 ^ for z e D, 

where L 1 ^ denotes the (generalized) Laguerre polynomial of degree q — 1 with 
parameter 1 . It is well-known that the Laguerre polynomial L^ j has exactly 9 — 1 
strictly positive roots, which implies that the Berezin density will exhibit a typical 
Fresnel-type ring pattern. This resembles what happens in the one-dimensional 
GUE case, where the zero density points for the Berezin density come from the 
zeros of the sine kernel. Actually, the analogy is more than a superficial similarity. 
If we consider rather big values of q, and scale down to local distance (mq)~^ 2 , 
with 

B^S') = (mq)- 1 B^(z + (mqy 1 ' 2 ?), 
then by the above we have, for zeD, 

fi&,<9 — q-'LUq-Wfe-^ 2 , 



asm,tn +00 with n = m + O(l). Next, if we let q — > +00, we get that 

y (-lmf' f = 7i(2|«'|) 2 
to f!(f + 1)! ' 
where ]\ is the standard Bessel function. The identity 



q- 2 l}Aq- x m z fe 



'12 



f. 

Jd 



,2iRe[{'Q 



cL4(Q 
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shows that we are dealing with a planar analogue of the sine kernel (the sine 
kernel is the Fourier transform of the characteristic function of the interval [-1, 1], 
the one-dimensional unit ball). 

We also investigate the local behaviour of the Berezin transform when |z| = 1, 
i.e., when z is on the boundary of the bulk. Using the same blow-up scale as with 
an interior point, we show that the blow-up Berezin density f>^„(£) tends to 
a limit which can be expressed in terms of Hermite polynomials (see Theorem 
5.10) . Here, too, there is a ring-like pattern in the interior direction, but it is not so 
pronounced as it is for interior points (the bulk). We express the 1-point intensity 
near a boundary point in terms of a sum of squares of Hermite polynomials. The 
Wigner semi-circle law then gives the asymptotic behavior of the 1-point function, 
which tells us the intensity of the process. We find that for big q, but much bigger 
m, n with m = n + O(l), the 1-point function is nontrivial in the annulus 

1 - q 1/2 m- 1/2 < |z| < 1 + q V2 m- 1/z ; 

inside it is essentially constant, and outside it more or less vanishes. Near the 
outward boundary of the annulus at the scale (mq)' 1 ^ 1 , we expect the statistics of 
the point process to be related with the well-known Airy process. 

Lifting to two complex variables. Analogous results to [2] are obtained on com- 
plex manifolds by Berman 1 6 1 . We note that the polyanalytic Ginibre processes can 
also be viewed as processes in C with the rather singular weig hte- |Zl|2 <5 (zi-z 2 ), 
where So is unit point mass at 0. Berman considers reproducing kernels of polyno- 
mial subspaces as the total degree of the polynomials tends to infinity. In contrast, 
here we discuss the case where one variable has bounded degree and the degree 
of the other variable goes to infinity. 

The polyanalytic Ginibre ensemble and weighted interpolation. It is well known 
in the theory of random matrices that 

1 1 

— det[K m/)U (z,,z ; )]" = \A(z lr ...,z n )\ 2 , 

where A(zi,...,Z„) = n/i : j<;(Zj - z,) is the van der Monde determinant. A point 
configuration in a compact set which maximizes the van der Monde determinant 
is known to be a good choice of nodes for Lagrange interpolation [21]. Instead of 
considering points confined to a compact set, one can add a weight to prevent the 
points from going off to infinity. This leads to the same expression which arises in 
the context of random eigenvalues. 

We turn to the polyanalytic case. One shows that 

1 1 

— — det[K m , M (z; ; z ; -)]"? j = lA^zi,...^)! 2 , 

where Z,„ „ is a normalization constant 

Z m ,n A := f |A 9 (zi, ...,z„,)| 2 e- mZ M Q(z ' ) cL4(z 1 ) . . . dA(z^) 

JC"i 
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and 



^(zx,...,z m ) = det 



( 1 
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Zl 




yil—l 
^1 


~tt— 1 






ZiZi 


ZnqZnq 


_fl-l 

z i 


7 


=1-1 n-1 
V ^1 z l 





So, A ? (Zi, ...,Z ncj ) is the polyanalytic analogue of the van der Monde determinant. 
As in the case q = 1 (which gives the usual van der Monde determinant), the 
expression |A 9 (zi, z„ ? )| 2 measures how good the configuration is for Lagrange 
interpolation by polyanalytic polynomials. So, the polyanalytic Ginibre ensemble 
is a way to produce random Lagrange interpolation sets, using the Gaussian 
weight for confinement. 

Further results and open problems. In Il20l and (3J, the authors showed that the 
fluctuations of eigenvalues of random normal matrices tend to Gaussian free field. 
The fluctuations of the polyanalytic Ginibre process will be discussed in a separate 
paper - the limit is again the Gaussian free field, but the variance depends on the 
degree of polyanalyticity. 

One could naturally address all the questions discussed in this article with 
a more general weight Q. We conjecture that the spectral droplet will be the same 
as in the analytic case. It is also likely that the blow-up of the Berezin density 
at a bulk point z will have a universal limit, which we here computed to be 

2. Polyanalytic Bargmann-Fock spaces 

An orthogonal basis. We will consider the Bargmann-Fock space A 2 (C) of poly- 
analytic functions of degree < q — 1, i.e., functions of the form 

9-1 

f(z) = JV/ r (z), 

r=0 

where all the components f r are entire, subject to the norm integrability condition 

H/H V) := H/I^Ce-^) - l \m\ 2 e- ml42 dA( Z ) < + oo. 

Here, as always, m > 0. We note that A^C) is the standard Bargmann-Fock 
space. As before, let Pol„ !ilJ/i; be the closed subspace of A 2 i(? (C), defined by the 
condition that all the components f r are polynomials of degree less or equal to 
n-\. Moreover, let X„ !/? and K mi „ rCI denote the reproducing kernels for A*, , ; (C) and 
Pol,„ „ /(? , respectively. We will be concerned with the asymptotic behaviour of the 
kernel K m ,n,q, as m, n —> +oo while n = m + O(l). 
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Bargmann-Fock-spaces of polyanalytic functions have been considered in, 
e.g., [25], where the reproducing kernels and orthonormal bases were identified. 
We will attempt to supply a self-contained account of these basic matters. 

To begin with, we identify an orthonormal basis for the space Pol m<H/(? . Here, 
and later as well, we will need some standard properties of the classical orthogonal 
polynomials. For details, we refer the reader to |16J. We need the generalized 
Laguerre polynomials ((x)j := x{x + 1) • • • (x + j — 1) is the Pochhammer symbol) 

Proposition 2.1. For q < n, the following functions form an orthonormal basis for 
Pol m/B< q (i, r, j, k are all integer parameters): 

e) r {z) := m (i+1)/ V V r (m\z\ 2 ), 0<i<n-r-l, Q<r<q-1, 

e 2 k (z) := ^ m (k+l)l2 z k L k (m\z\ 2 ), 0<j<q-k-l, l<k<q-l. 

Proof. Clearly, all the above functions belong to the space Pol m/fl/? . Also, all the 
functions ej r are orthogonal to each of the functions e 2 k , for the indicated ranges 
of the indices, as can be seen by integrating along circles using polar coordinates. 
Next, we show that the functions e. form an orthonormal set. Any two functions 
there having different parameter i are orthogonal, again by integrating along 
circles using polar coordinates. So, we fix i, and pick two indices T\,r%. We compute 
the inner product of two such functions: 

f e i ) . i (2)^)e- m ' z ' 2 cM(z) 

= 1 l z l 2, ^(-l z l 2 )^ 2 (-l z l 2 )e- m ' z|2 dA(z) 

where the delta is in Kronecker's sense. In a similar fashion, the functions e 2 , 
form an orthonormal set. So, the functions e 1 ,e 2 , together form an orthonormal 

' ;,r ],k o 

set. Next, the dimension of the span equals the total number of vectors, which we 
calculate to nq, which equals the known dimension of the space Pol m „ The proof 
is complete. □ 

The reproducing kernel. For q < n, we conclude that the reproducing kernel of 
Yo\ m ,n, q equals 

n-r-l j 

(2.1) K m ^(z,w) = mY J Yj 7^(mzw)'Ll(m\z\ 2 )V r (m\w\ 2 ) 

r=0 1=0 ^ '' 



9-2 9-Jr 1 ■, 



;'=0 k- 
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We note that by plugging in n = +00 in Proposition 2.1 we get an orthonormal 
basis for the space A^„(C). It follows that the same procedure of plugging in 
n = +00 in the above expression for K nii , h q gives us K„, „, the reproducing kernel 
for A 2 ncj (C). What is probably less obvious is that K Mi q may be written in a much 
simpler form (this representation is, however, known; compare with [5], lfl9l ). 

Proposition 2.2. We have that 

<?-l oo ( 

K m/Ci (z,w) = m^Y^ , r + a, (mzwyV r (m\z\ 2 )L l r (m\w\ 2 ) 



r=0 i=0 
9-2 '7-7-1 



+ m Yj £ ^^(mz^L^mlzftL^mN 2 ) = m L^m^ - w\ 2 )e" 



j=0 h 



Proof. It is clear from Proposition |2.1| that the double sum expression equals 
K mj q(z, w), so the only thing that needs attention is the last equality. We first do the 
case w = 0. Many of the terms in the double sum vanish, so that we are left with 

q-l 

X, M (z,0)=m^L°(m|z| 2 )L0(0), 

It is well-known that L°(0) = 1 and that JL'lJ) — so that the above reduces 
to 

K, M (z,0) =mL]_ l {m\z\ 2 ). 
We turn to general w eC.lt is well-known that the transformation 

T w [f](z) := e -i>«M 2 ->«^/( z + W ) 

acts unitarily on A 2 ^(C); its adjoint is T*„ = T- w . By the reproducing property of 
K m , q (z, 0), we have, for / e A 2 (C), 



-1«,|7:,|2 



e - 2 «nm f{w)=Ta[f]{0)= Ta; /(z)K ffl , £ ,(z,0)e- m|2| dA(z) 



-l 

= f /(z)TUX, M (-,0)](z)e- m ' z ' 2 dA(z) 
Jc 



= OTe -|»N 2 /( z )Li(m|z-zi;| 2 )e fflZll 'e- ffl|2|2 cL4(z), 
Jc 

The claim that K miCj (z, w) =m Lj_ 1 (m|z - w\ 2 )e mzw now follows immediately. □ 

Szego asymptotics. We shall show that the kernel K m ^ /Cj {z,w) is approximated 
well by K nii q{z, w), provided that z, w e D are rather close to one another. It will be 
instrumental to study the asymptotic behavior of the partial sums 

Em = t jfi 

of the Taylor expansion of the exponential function. In [23], Szego showed that 
£„(«) = , _ 1 (fel - t)t C + , 



(2 - 2) - - VSJE- M-C 
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provided that ICe 1 e | < 1 and \Q < 1, whereas 

(23) = 1 (Ce^)^(l + e£(Q), 

provided that ICe 1- ^ < 1 and |CI > 1. Here, we have the convergence 

lim el(z) = lim e?(z) = 

uniformly on compact subsets of the respective domains. As for \2.3\ , the uniform 
convergence e? (Q — > as A: — » +oo holds also in certain unbounded subdomains; 
in particular, along the real line, we have uniform convergence on all intervals 
[a, +oo[ with a > 1. 

An elementary estimate of Laguerre polynomials. The following elementary 
estimate of generalized Laguerre polynomials will prove useful. 

Lemma 2.3. Suppose a is a positive real. Then, for k = 0, 1,2, . . ., we have that 
\L a k {x)\ <l~{x + a + kf, x e [0, +oo[. 

Moreover, with f, := a + 2k - 2 + 2 ^\+{k- l)(k + a - 1), zt>e foawe f/zaf 

~(* - /8) fc < (-l)^(x) < I**, x e \0, +oo[. 
Actually, the related inequality \L?(x)\ < jyx* holds for all x 6 [|jS, +oo[. 
Proof. We begin with the estimate (0 < z < k is assumed) 

(a + i + l)fc_; < (a + fc) ; . 
Next, we note that for x > and = 0, 1, 2, . . ., we have 

> Z-i v > i\(k-i)\ L-> i\(k-i)\ Lu. i\(k-i)\ 



m-i)\ Lu iHk-i)\ i-i il(k-i) 

i even I odd 

and in view of the above estimate, 

i even i even 

and, analogously, 

E ( ^^Et^ = ^(<-*«>'-<-*-*>')- 

/odd v ' i odd v ' 

By discarding alternatively the even or odd contribution, we arrive at 
\L«(x)\ < ^-((x + a + k) k + \a + k- x\ k ), 



2kV 

which is slightly stronger than the first estimate. 

It is well-known that L£ is a polynomial of degree k all of whose zeros are 
distinct and real, and they all fall in the interval ]0,jS[ (see 18.16.13 of [16|). The 
second estimate follows immediately from this. □ 

Approximation of the polynomial polyanalytic reproducing kernel. We estimate 
the difference K m , M - K„, /q . 
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Proposition 2.4. Let z,w e C be such that \zw\ < B\ < 1 and \zw\^ |ziy| < £?2 < 1- Lef 

Mbea positive real number. Then, as m, n — > +oo ztfe'Ze \m - n\ < M, 

\K m , n ,„{z,w) - K m , q {z,w)\e- m ^ < Cm"e^(l - d^l + \z\ 2 * + \w\ 2 \ 
where the constant C depends on q and M. 

Proof. In view of ^2.1) and Proposition [Z2j we have that for q <n, 

<?-l +oo i 

K mM {z,w)-K m>nM {z,w)-m'Y_ l E t- ' -x, (wzi(;)'L;(m|z| 2 )L^(m|a;| 2 ). 

r=0 i=H-r ' 



By Lemma 2.3 we get that 



(2.4) 



E 7^(^sO'^|z| 2 )LXmM 2 ) 
(r + z)! 

< ^ — — (m\zw\)'(m\z\ 2 + z + rY(m\zv\ 2 + i + r) r 

i=n—r 
1 + °° 1 

= — 7 — (m\zw\) n ~ r+1 (m\z\ 2 + n + j) r (m\w\ 2 + n+ if 

r\ *—> (n + /) 

< — E (—^(mlzMr-^mzfy + (n + jf][(m\zv\ 2 Y + (n + jf]. 

For r confined to0<r<^-l and for big n, say n > no(q), we have 

(n + j) 2 * 2 (n+jf 2 
< ■ — — — and — < 



(n + /)! (n + j- 2r)\ (n + ;')! (n + j - r)\ ' 

so that if we use the notation 



Rk(z) = £ \z> =e z - E fc (z), 
j=k+l 1' 



we obtain from | |2.4| that for < r < q - 1 and n > Ho(q), 



(2.5) 



r! 



E T^A™w)%{m\z\ 2 )V r {m\w\ 2 ) 
*—i (r + 1)\ 



(r + i) 

—i 

4T-i m 



i=n—r 
iX-l m r 

< 



f [|zwrR»_i(m|zw|) + 2(|z| 2r + \w\ 2r )R n - r -i(m\zw\) + 2|zzt>| r R„_ 2) -i(?w|zw|)]. 

Next, we see from Szego's asymptotical expansion \2.2\ that for large k, I with 
I = k + O(l), we have 

(2.6) = (2^ W^Cze^Ya + <£<«:/*)), 

as k, I — » +oo and / = A:+0(1). By a careful application of | |2.6} to \2.5\ , and summing 
over < r < q - 1, the assertion of the proposition follows. □ 

Corollary 2.5. For z,w 6 D, fef 1 - |zzc| > t > 0. TTzen, as m,n — > +oo, wMe 
|m - n| < M, 

K^few) = K^frw) + 0(e-5"' T2 e m|zi " l ) / 
zfftere f/ze "O" constant depends on x, q, and M. 
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Proof. A Taylor series expansion of the logarithm gives that 

te 1 -* <e-$V- t)2 , 0<f<l, 

and, together with the fact that the exponential function grows faster than any 
given power, the assertion follows from Proposition |Z4] □ 

As we shall see, Proposition |2 .4| implies that the Berezin density 



\K mr n, q (z,ZV)\ 2 



n,q 

behaves locally near z like the Berezin density 

\K m Jz,w)\ 2 , |2 m 2 L\_ x {m\z-w\ 2 ) 2 e 2m ^ z *> 
( ] ^ (W) ~ K m , q (z,z) ~ ml* (0)e»W 

= -L\ Jm\z-w\ 2 ) 2 e- m]z - m]2 
q i 

To make this precise, we recall the definition of the blow-up Berezin density: 
Proposition 2.6. Fix z e D. Then 

|bL%(5) - fX-i (l£l 2 )V ia2 M(£) — » 0' 



asm,n — > +oo zp/zz/e n = m + O(l). 
Proof. By l [2.7) , we have 



fi«(0 = m^E^z + m" 1 ^) = <? - 1 Lj_ 1 (|5| 2 )V^, 

so the comparison is with the blow-up Berezin density for K. m n. We write w := 
z + m~ l l 2 li; since z e Dis fixed, we have 1 - \z\ 2 > z > for some small t, and we 
suppose that w is close to z so that 1 - \zw\ > t > as well. Then, by Proposition 
2.2 and CoroHary |2.5{ 



\K„,^(z,w)\ 2 _ mN2 = \K miq (z,w)+0(e-2™e m W)\ 2 
K m , M (z,z) 6 " K m/ ,(z,z) + 0(e-^ 2 e <*l 2 ) 

2 



ImL^^mlz - w\ 2 )e mzw + 0(e-i mT e m|ZB ")[ 



(mq + 0(e—2 nn ))e"'W 

so that if we we expand the square using that 

\L]_ l (m\z-w\ 2 )\ = 0{m^), 

where the "O" depends only on q (this follows from Lemma [23) , and simplify the 
expression, we arrive at 

B< z> „» = -L 1 Am\z- w\ 2 ) 2 e- m]z - w]1 + 0(m^V^e^«- |, * l ) / 
' '' q 1 

which immediate gives that 

b£^(w) = -L 1 ^(nilz - w\ 2 ) 2 e-' nlz - wl2 +0(wj 2< ?- 2 e -5 mT2 ). 
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The corresponding blow-up Berezin density then has 

(2.8) fi^tf) = m-^%,{z + m-^) = ^L^fe^ 2 + O^e"^). 

As exponentials grow faster than polynomials, the error terms is negligible for 
big m. For fixed zeD, the requirement on £, so that 1 - \zw\ > x for some fixed 
t > is fulfilled for big m if, say, |£| < log m is required. So, | |2,8) has the immediate 
consequence that 

(2.9) f \£%„/Z) - q-V J\Z\ 2 ) 2 e-W\ dA(£) = oUMlog/^V^" 2 ), 

JD(0,logm) 

where more generally D(zo, p) denotes the open disk of radius p about Zq. Since 
the associated blow-up Berezin measures 

dB^) := fi^( W )dA(£), dB« (€) == ^L^ftV^dAtf), 

are both probability measures, the assertion of the proposition follows from | |2.9) 
once it is noted that the right-hand side of < |2.9) tends to as m — » +oo. □ 

We note that the claimed convergence dB^L — > dS z forz e D is an immediate 



consequence of Proposition 2.6 



3. Berezin density asymptotics for an exterior point 

Convergence to harmonic measure. We show that the Berezin measures have the 
convergence dB^„ (; —> da> z for zeUasm^-* +oo while n = m + O(l). Here, co z 
is harmonic measure with respect to the point z and the exterior disk D e . 

Concentration of the Berezin mass. We first study the concentration of the Berezin 
measure to neighborhoods of the closed unit disk D. We recall the standard 
notation D(zo, p) for the open disk of radius p centered at Zq. 

Lemma 3.1. Suppose z e D c and that p > 1. Then 



(3.1) f E% n/q (w)dA(w)^0, 

Jc\D(0,p) 

as tn, n — > +oo w;fa7e n <m + O(l). 
Proof. Let w e C \ D(0, p). We note that 

(3.2) n - 1 < m + O(l) < 2m 



:<z> 



for sufficiently big m, and so, by Lemma 2.3 

n—r—l . 

< 3 ' 3 ) E r^\^ mz ^(Mz\ 2 )V r (m\w\ 2 )\ 
(r + 1)1' 1 



n-r-1. ^ 

< 7 — — (m|zwiy(m|z| 2 + i + r) r (m\w\ 2 + i + r) r 

L-i rUr + lY 



rl(r + i)\ 

!=0 V ' 



- | % 7^(m\zw\Y +2r < 9 -(m\zw\Y E^imlzzvi: 



(r + i) 

(=0 v ; 
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If we plug in n = q in 1331, we obtain 

(3.4) J_j -2—\{mzwfL)(m\z\ 2 )L)(m\w\ 2 ^ < -(m^E^m^). 

k=l V + '' 1' 

We may restrict to q < n; after, we are considering the limit as n — » +oo. As 
E 9 _i < £„-i on the positive half-axis for q < n, an application of ( |3. 3| > and < |3.4) to 
( |2l) gives 

(3.5) |K„ I/M (z,a;)| < 2me 9 {m\zw\f- l E n -i{m\zw\). 

Next, we observe that Pol m , nj i C Pol m/ „ # (since q > 1), which implies that 

K OT , M (z,z) > K m , n ,i( z ' z ) = mE n -i(m\z\ 2 ). 
We conclude that the Berezin density may be estimated as follows: 

(3.6) Bt) n ,{w) < ^\m\z W \f^ E ;- l{ f Z ^ e 

' E„-i(m\z\ z ) 

Finally, we see from Szego's asymptotical expansion \2.3\ that 



m\w\- 



(3.7) = (inkr^Wf^itf-Wni + 4(K/k)). 

k < I + O(l), where the convergence e 2 (lC,/k) — > is uniform if C is real with C, > a 
for some fixed a > 1 . This leads to 

= ^ n - D] _1/2 K» - i)Anrv-vr 2 a + o(D), 

E„_i(m|z| 2 ) 

where the "o" term is uniform in the convergence. As we implement this estimate 
in | |3.6| , and integrate over C \ D(0, p), the claim follows. □ 

A principal value calculation. We follow the approach of |2], and calculate a 
certain principal value integral. 

Lemma 3.2. Fix z e D e . For any I = 0, 1, 2, . . ., we have 

pv I i^E^(w)dA(ie) -» z-\ 
Jc 

asm,n — > +oo jcif/i n = m + O(l). 

Proof. The case ^ = 1 was treated in J2J, so we may from now on assume that 
q > 2. We write 

Y — fe"' 4- K"^ 

<^m,n,q — ^m,n,q ~ *^m,n,q> 

where 

1-1 n-r-1 . 

< /M (z, W ):=m£ £ -J^(rnzw)%{m\z\ 2 )K(m\w\ 2 ) 

and 

,-2 q-y-l 

f^ M (z,o;):=m£ £ -^{mzw) k L){m\z\ 2 )L){m\w\ 2 ). 

j=0 k=0 " 

It follows that the expression \K m A /? | 2 decomposes accordingly: 



(3.8) \K m , n , q (z,w)\ 2 = |X^(z,^)| 2 + |< M (z,^)| 2 + 2Re[K^(z,a;X M (z,H;)]. 
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We first consider the following integral involving \K ! mn(j \ 2 : 

(3.9) pv f w-<\K^ n /z,w)\ 2 e- m W 2 dA(w) 
Jc 

q-l n-n-1 h-i'2-1 „ | | 

E> 7 pv I w~ l , ' — —(ntziv) h (mzw) 
f^n ft! (r! + ii)!(r 2 + z 2 )! 



= m 

rij2=0 h=0 i2=0 

xLX{m\z\ 2 )L\{m\w\ 2 )L^{m\z\ 2 )^ 

q-l n-n-ln-r 2 -l ; 1+ , 2 . . 

- V 'E E E 

ri,r z =0 h=0 i 2 =0 v 1 ly ' v z z; ' 



x IwrL^wilwl^L^mlrol^e-"' 1 ^ dA(a>), 
Jc 



fl 

where the delta is understood in Kronecker's sense. The identity (for p = 1, 2, 3, . . .) 
(3.10) Lr P (x) = g( r "^ P 1 " 1 )L«(x) 

plus the standard orthogonality properties of the Laguerre polynomials gives that 



r*+oa 

Jo 



„ m( fHt _ (ri+ii)!(r 2 -r 1 + Z-l)! 

o r a !(r 2 -ri)!(Z-l)! ' 



where the right-hand side should be interpreted as for r 2 <T\. By polar coordi- 
nates, then, we have 



\ \w\ Zh L' r \(m\w\ z )Ll\ + \m\w\ z ) e - mwr dA(w)=m 
Jc 

and < |3.9[ simplifies to 

(3.11) pv f w-'\K' m (z,w)\ 2 e - m]l " ]1 dA(w) 
Jc ' 

n=0r 2 =ri i 1= o v ; v ; v ; 

provided n is so big that q + l <n. Next, we apply Lemma 2.3 and ( |3.2| > (using that 
z e D e ) to arrive at 

n-/-r 2 -l . . 



< 3 - 12 ) E (r rfeni (^lzl 2 )' 1+ '|^(^lzl 2 )4 +i (^lzl 2 )| 
(r 2 + n + /)! 

n-l-rz-l 

* E ? I • I m («|z| 2 ) il+ '(«|z| 2 + »i + n) n (m\z\ 2 +h+l + r 2 ) r 

r—i \1"2 + h + ')■ 



!'l=0 

n—Z— 72—1 



^ y i+ri Y t — 5 — ™ Nzi 2 r +, ' 2+,i+; < 3'' i+, ' 2 c«izi 2 )'- i £„_i(ff2izi 2 ). 

(r 2 + 
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On the other hand, by the estimate from below in Lemma 2.3 

n-r-l 

(r + i)! 



(3.13) K m/M (z,z) > < <M (z,z) = m£ £ -^ {m \z\ 2 )%(m\z\ 2 ) 

r=0 i=0 



r\(r + i)\ 

r=0 i=0 v ' 

provided that m\z\ 2 > f5(n), where 

/3(n) :=n + q-4 + 2 J\ + (q- 2)(n - 1) = n + 0( Vn). 
As we assume z e D'' and n - m + O(l), we must have 

ffi|z| 2 - j8(n) = m(|z| 2 - 1) + m - jg(n) > |m(|z| 2 - 1) 



for big enough m, n, and so, by < |3.13[ l, 

?-l n-r-l 

r\(r + i) 



(3.14) K m , M (z,z)>K^(z,z)>m£ £ W^^ 2 ^ 



r=0 i=Q 



-(mW-^En-ximW 2 ) - E q . 2 {m\z\ 2 )} > — {mW^E^m^ 2 ). 



~ (<7-l)! v v ' ' ' r y (9-1)1 

Here, we used that the E n -i(m\z\ 2 ) is much bigger than E,,_2(ffz|z| 2 ) as m,n both 
grow. Let us look at the contribution from < T\ < q - 2 in the right-hand side of 
( [3.11) using the estimate \3.12) : 



r 1= 0r 2 =ri i 1= v z l; v ; v z 1 ; 



ri=0i- 2 =r 1 v 17 

^md^O^lz^^E^^mlzl 2 ), 



for an appropriate positive constant C\{q,l). As we combine this estimate with 
( 3.14) , we obtain 



|z| y y y r 2 \(r 2 -r x + \ -\)\ 

K m , n , q (z,z) L> Lf U (r 2 - n)!(I - l)!(r 2 + ij + /)! 



ri=0r2=ri ^=0 

x (m|z| 2 ) ! ' 1+ '|4(m|z| 2 )4 + '(m|z| 2 )| < C 2 {q,l)m~ x = 0{m~ l ) — » 0, 

as m,n —> +00 in the prescribed fashion. So the contribution to \3.11\ which comes 
from < T\ < q - 2 is negligible from the point of view of the Berezin density. It 
remains to consider the contribution from r\ = q - 1. The corresponding part of 
the sum in the right-hand side of < |3. 1 1| > equals 



n-l-q 



, =0 fo + ii + Z-1) 
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and we now claim that 



(3.15) m 



z ! "P (<?-!)! 



(a + u +. 



K m , n Jz,z) 4—* (q + h + 1-1)1 

/}— U 

as tn, n — > +oo in the given fashion. We use the recurrence relation [ 3.10) to write 



9-2 



(3.16) 



L l ;_ 1 (x) = L i £(x)-Y J 

s=0 



q+l-s-2 
1-1 



Q(x). 



A straightforward argument allows us to show that as insert this into 13.15) , the 
sum that is subtracted on the right-hand side o f (|3.16) makes asymptotically no 
contribution to the sum in p. 15) . Consequently, p. 15) is equivalent to having 



n—l—q 



(3.17) 



m 



(?-!)! 
(q + h +1-1)1 



as m, n — > +oo in the given fashion. If we insert the expression f2.1) defining 
K m , n ,q{z,z), we see that (there is some cancellation of terms) 



n-l-q 



(3.18) ^-m^^ £ (? ;?; + /l 1) , (m^^>N 2 )^>^) 



<j-2 



;w E E t^N 2 ! 2 )^!*! 2 ) 2 

,=0 i=0 ^ '' 

4=?tH ft £4 (o-l)' 



K mi n r q{z, z) 

+ m L E (jTXji^l 2 ! 2 )^! 2 ! 2 ) 2 + ™E (g +^-i)! < OT ' z ' 2 )'' L i-i< OT ' z ' 2 ) : 



ii=0 



By careful application of Lemma 2.3 to all the Laguerre polynomials in this ex- 
pression, while inserting 13.13 1 to control the denominator, we indeed get 1 3.17 1. 
So, after a lot of effort, we have obtained that 

(3.19) pv w~ l ' — — e""' N &A{w) — -> z~ l 

as m, n —> +oo with n = m + O(l). Analogous but slightly easier arguments (left to 
the interested reader) show that 



K„/Z,TV)\ 2 



pv w 



dA(w) 







and 



^L,n,n( z ' w Wm,n,q{ z > w ) 



e~^dA{zv) — > 0, 



^m,n,q( z / z ) 

again as m, n — > +oo with n = m + O(l). Finally, we put everything together based 
on the decomposition ||3.8): 



\K milhq {z,w)\ 2 



-\wf 



JC JC K- mn q(z,z) 

as m, n — > +oo with n = m + O(l). This ends the proof. 



dA(o>) 
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Convergence to harmonic measure. We now show that even in the principal 
value sense, the Berezin density tends to avoid the interior of the unit disk. 

Lemma 3.3. Fix a real parameter p with < p < 1 and a point z e D e . Then, for 
I = 0, 1, 2, . . ., we have the convergence 



pv f ut'eJ® (w)dA(w) — ► 0, 



'B(0,p) 

asm,n — > +oo zpjf/j n = m + O(l). 

Proof. In terms of the decomposition ( |3.8L we will focus on the the term \K^ n „\ 2 
and leave the other two to the reader (the necessary arguments are similar but 
slightly easier). The analogue of |3.9) reads 

(3.20) pv f w~ l \K] n (z,wfe- m ^ 2 dA(w) 
Jo(o,p) 

(\— 1 min{n— ri— l,n— Z— f2— 1} ^ ^ 

=- v 'E E ~ (, 1+ J;^, + /). ^ |z|2 )' 1+;l ^ ot|z|2 ^ + '^ |z|2) 

n,r 2 =0 i 1= v 1 11 y 1 1 ' 



x f (m|o;| 2 )' 1 L^ i (m|w| 2 )4 + '(m|ii;| 2 )e- m|l ' ;|2 dA(a;). 

Jo(0,p) 



By Lemma 2.3 we have 

f (OT|w| 2 ) !l |L^(m|i(;| 2 )L^ + '(OT|it;| 2 )|e- mNJ dA(iy) 



< — *— f (ffzM^'^wzM 2 + fj + n) ri (m|af +i t + 1 + r 2 Y 2 e-' mwr dA(w) 
nW- J D (o,p) 

2»'i+''2-2 



+r 2 -2 

< - (mN / ) ,I [(OTN z ) ri +(fj +ri) ri ][(w2N z )'' 2 +(h +/ + r2r 2 ]e- m|ar cL4(a;) 

r l- r 2- Jo(0,p) 



2 ) il [(m|o;| 2 )'' 1 + (h + r 1 ) ri ][(m|w| 2 )'' 2 + (h ->- ' -^i -mM 2 , 

2fi+r 2 -2 ( rmp 2 nmp 2 



>-2 ( pmp- pmp- 

— \ f il+ri+r2 e- f df + (n + n) ri t ,1+,2 e- f dt 
!w I Jo Jo 



ri!^!??! 

Xmp 2 r*mp 2 \ 

In terms of the function 

1 r 6 

*(«,&):= r , fe-'dt, a,be[0,+oo[, 

r(fl + 1) Jo 

which takes values in the interval [0, 1 [, the estimate becomes 
f (mliopJ^lL^ml^ftL^CmM^le-^dACH;) 

JlD(0,p) 

2i-i+r 2 -2 r 

< r !r !m |( z 'i + ri + r 2)!AT(ii +n + r 2 ,mp 2 ) + (h + n) n {h + r 2 )lx(h + r 2 ,mp 2 ) 
+ (z'i + I + ^Hzi + riMz'i + n, mp 2 ) + (h + riHzj + / + r 2 ) r2 hlx(h, mp 2 ) \. 
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and if we use that a i-» xi a > b) is decreasing for fixed b (a direct calculation involving 
derivatives suffices to verify this), we get 



(3.21) f (m|tDp) fe |l4(m|tD| 2 )l4 +I (m|o'l 2 )|e~ w,ls|a dLA(H7) 

JD(0,p) 

2ri+r 2 -2 



(ij +n +r 2 )! + (fj +r 1 ) ri (/ 1 +r 2 )! 



+ (ii + / + r 2 )'' 2 (z'i + ri)! + fa + r^fo + I + r 2 ) r H^x(h,mp 2 ). 

Next, since z e D e , n - m + O(l), and i-y + 1 < n — 1, we may use another aspect of 
Lemma |231 to see that 



(3.22) \V r \(m\z\ 2 )L\} 2 + \m\z\ 2 )\ < -L( m | z f)^, 

provided m, n are big enough. As we combine the equality | |3.20) with the estimates 
( |3.21) > and < |3.22| , and use some well-understood comparisons of factorials and 
powers, we arrive at 



(3.23) 



pv f w- l \Kl (z,w)\ 2 e- m W 2 dA(w) 
Jd(0,p) 



n-l-l 



<mC 3 (c,,l)\z\-'^- 2 £ tpj— x{h , m p 2 ), 



2\h+l 



(h + /)! 

for some appropriate positive constant C^{q,l). The function x(h,ntp 2 ) drops off 
exponentially quickly to as i\ exceeds mp 2 by a margin greater than 0(m 1 ^ 2 ), as 
can be seen, e.g., by an application of the Central Limit Theorem (compare with the 
next section). This means that effectively we are summing up to m\p\ 2 + 0(m 1/ ' 2 ) 
in the right-hand side expression of | |3.23), which does not permit the sum to 
compare with the size of K„, /Hj(? (z,z); cf. ||3.14[|. This results in the convergence 



r \¥} mn {2.,wf 

P v vr l ^ . — r e dA(w) — > o 



as m,n — > +oo with n = m + O(l). Together with the estimates which were left as 
an exercise to the reader, we get 

pv r w _, \k {z w)\2 e _, m2dA{w) _^ o 

Jo(0,p) K m ,n,q{Z,Z) 

asm,n — > +oowithn = m + 0(l), which amounts to the assertion of the lemma. □ 

As in the proof of the Theorem 2.11 in [2 J, we may now conclude the follow- 
ing. 

Theorem 3.4. Fix z e D e and a bounded continuous function g on C. Then 

as m,n — > +oo with n = m + O(l). Here, dco{w,z,W)*) is harmonic measure with respect 
to the point z and the domain D e . 
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4. Poly-Bargmann transforms 

Purpose of the section. In this section, we discuss the poly-Bargmann transforms, 
a generalization of the classical Bargmann transform, which are needed later 
when we analyze the Berezin density at a boundary point. The poly-Bargmann 
transforms appeared in Vasilevski's paper l25l , where the basic properties were 
presented. 

The Hermite polynomials and the Bargmann transform. We denote by Hi the 
j-th Hermite polynomial with respect to the Gaussian weight e - ^ 2 ("probabilistic 
Hermite polynomials"). The generating function identity 



7=0 

allows us to write 



1 c zMzH' 2 = y 2A 

^) 1/4 fef (271 



Hj(t)e-^ z i 
x 



( 2 ^) 1/4 j^W4f- # 

We recall the standard definition of the Bargmann transform: 

B[/](z) = f e zl -^ t2 f(t) dt, f e L 2 (R). 

As the function systems 

^— pHy^e"** 2 and 

form orthonormal bases for L 2 (R) and the Bargmann-Fock space A? -(C) (this is 
A^(C) with m = q = 1), respectively, we obtain the following well-known fact. 

Proposition 4.1. Tfe Bargmann transform B : L 2 (K) — > A 2 j(C) acfs isometrically and 

bijectively, and for each j - 0,1,2,..., the basis function (2n)~ 1/4 (j\)~ 1/2 Hj(t)e~i t2 is 
mapped to the basis function (j\)~ 1/2 zL 

A class of auxiliary operators. Let d z , d z denote the standard Wirtinger differential 
operators 

d z := \{d x - id y ), d z := l(d x + id y ), where z = x + iy. 
For r — 0, 1,2, . . ., we introduce the operator 

T r [/](z) := ^e' z ' 2 ^j/(z)e-' 2 ' 2 }, 

with the semi-group property Ti o T r _i = r 1/2 T r . We also consider the dilated 
variant 

W/](z) = T r [f m -v2](m^ 2 z), 
where f m -m(z) = f(m~ 1/2 z). It has the semi-group property T,„ ; i o T m , r _i = r 1/2 T mjr , 
and may be expressed in more concrete terms: 

W2 r ^ 2 mW\1 -a, I v -m| z | 2 



W/](z) = — e m ' z ' 2 ^j/(z)e 
We now study the effect of T,. on the basis elements (;'!) _1/2 zA 
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Proposition 4.2. For j > r, we have 



r z> 1 






= 4^ 



while for j < r, 



z-' 



Vr! 



ze^VJe-W 



z-' 



,}" 



(-1) 



^'-'Lnizf) 



Proof. The proof is based on an induction argument. The statement is obviously 
true for r = and all j. So, by induction, we assume that the statement holds for 
some r — 1 > and all j. In case j > r, we then have 



^{e- |z| V} = d z {(r - l)!z^ +1 L£ +1 (| Z | 2 )e- |; 

= (r - l)!z^{(; - r + l)L;i; +1 (|z| 2 ) - i£ +2 (|z| 2 )|z| 2 - |z| 2 Ln +1 (|z| 2 )}e-' z ' 2 

= r!Lr r (|z| 2 )z^''e- |z|2 , 

if we use the standard identity rL"(x) = (a + 1 - x)L"*}{x) - xL"* 2 (x). Next, in case 
j < r — 1, we have instead 



<%e- |z| V} = (-ir 1 -^{y!z r - 1 -a^ 1_i (|z| 2 )e- 1 

= (-l)'-^-!z^"jL;:((|z| 2 ) + Lp- 1 (|z| 2 )je- |z|2 = r-'L^lzlV^ 
which completes the proof. □ 

Pure poly-analytic Fock spaces. Write e ; (z) := (j\)~ 1/2 z> for / = 0, 1,2, . . ., which 
functions form the standard orthonormal basis for the space A 2 1 (C). The function 
T r [gy] e A 2 r+1 (C) is then a polynomial in z,z, where the degree in z remains 
equal to j, and the degree in z equals r. For general positive m, the functions 
et m {z) := (j\)~ 1/2 m^i +is> z>, with / = 0,1,2,..., form an orthonormal basis for the 
space A 2 n j(C). The functions T m A e j,m] are computed using Proposition 4.2 above, 
and we then recognize that we can identify them with the basis elements which 
appear in Proposition 2.1 We clearly have that 

span{T„ v [e //ffl ] : < /' < n - 1, 0<r<q-lj- Pol,„, M , 
and in view of the orthogonality properties in Proposition 2.1 we also must have 

(4.1) span{T m , r [ey <m ] : < /' < n - 1} = Pol„ w+ i e Pol w , 

where the "Q" is with respect to the inner product in L 2 (C, e~ m ' z ' 2 ). Following 
Vasilevski [25 j, then, we define the pure polyanalytic Fock space of level r + 1: 

(4.2) <5Pol m ,„,,. + i := Pol wl 9 Pol,„,„,,., 5A 2 m n r+1 (C) := A 2 >r+1 (C) 9 A 2 , r (C), 

with the understanding that Pol m/ „ o := {0} and A 2 Q (C) := {0}. The operator T m l - 
now becomes an isometric isomorphism 

Tm, r ■ A 2 j(C) -» bA 2 m n r+1 (C), T m , r : Pol m/ „,i -> 6Pol m ,„, r+ i- 
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We obtain the orthogonal decompositions 

q-l q-l 

Pol m , M = SPoUi, Ai/C) = 6A 2 mr+1 (C). 

r=0 r=0 

As a consequence, if K 6:mi „ >r+1 and JQ> :m , r +i denote the reproducing kernels for the 
spaces 6Pol m , n/r+ i and <5A 2 r+1 (C), respectively, we must have that 

q-l q-1 

(4.3) K m n q = ^ * Kd:m,n,r+lr ^-m,q = ^ ' t t^5:m,r- 

r=0 r=0 



The poly-Bargmann transforms. For m = 1, the poly-Bargmann transform of level 
r is defined by 

B r [f] := T r o B[f]. 

Proposition 4.3. We have 

B r [/](z) = ±— f /(t)H r (t - z - z) e tz -^ t2 dt, r = 0,1,2 

(2ti)1/4 J K 

Proof. We proceed by induction. For r = 0, the formula reduces to the usual 
Bargmann transform (if we recall that H = 1). By the induction hypothesis, we 
suppose therefore that the formula is valid for all integers up to r - 1. From the 
semi-group property Ti o T r _i = r 1/2 T r , we see that 

B r [f] = T r o Blf] = r 1/2 Tj o T r _j o B[/](z) = r 1/2 T a o B r _ x [/]. 

Since the formula holds for r - 1, the following calculation shows that it holds for 
r as well: 

r- 1/2 e |2|2 <9 2 ([(r - l)!]" 1 ' 2 ^* - z - z) e'H^V 212 ) 

= (r!)" 1/2 (-^(t -z-z) + (t-z- z)H r -i(t - z - z)) e^-i' 2 

= (r!)- 1/2 H r (f-z-2)e fe -5 z2 -J f2 ; 

here, we used the standard identity H r (x) = xH r _i(x) - fT^x). The proof is com- 
plete. □ 



5. Reproducing kernel and Berezin density asymptotics for a 
boundary point 

Purpose of the section. In this section, we will calculate the limit of the blow- 
up berezin transform at a boundary point z, that is, |z| = 1. There is no 
loss of generality to take z = 1. Our strategy is to investigate the blow-up of the 
reproducing kernel of the space of analytic polynomials Pol m , n/ i (with q = 1) first, 
and then use this information together with poly-Bargmann transform to lift the 
the asymptotics to the context of the general polyanalytic spaces Pol m ,„ /I? . 



24 



Haimi and Hedenmalm 



The central limit theorem revisited. The following improvement of the central 
limit theorem will be needed. Let cdfx denote the cumulative distribution function 
of a real-valued random variable X and let 

erf(z) = — L= f e'^df 
y2n J-oo 

be the error function. We shall write i.i.d. as shorthand for independent identically 
distributed in the context of random variables. The following result is from Q, 

eh 

Theorem 5.1 (Berry-Esseen). Let Xi,X2, ... be i.i.d. real-valued random variables with 
E(Xj) = 0, £(X 2 ) = 1 and £(|X ; f) = p < +co for all j. Also, let Y n = n _1/2 £<LiX/. 
Then there exists an absolute constant C such that 

IcdfWx) - erf(x)| < -p, x £ R. 
yn 

The Berry-Esseen theorem gives the following asymptotics for the partial 
Taylor sums of the exponential function. 



Lemma 5.2. We have 

E„-i(x) 



erff^U 0(rT 1/2 ), 
V vn / 



flsn-> +oo, uniformly in x e [0, +oo[. 

Proof. Let Xi,...,X„ be independent exponentially distributed random variables 
on [0, +oo[ with density e~*. It is well known that the sum ]£? =1 Xy obeys a gamma 
distribution with the cumulative distribution function 1 - e _;,: E„_i(x). The random 
variables X\ - 1, X„ - 1 all have zero mean and variance 1, and the third moment 
is finite, so by the Berry-Esseen theorem, we have 

(5.1) 1 - e-*E n _ a (x) = cdf z;UX/ (x) = cdf £ » =1 x r „(^) = erf(^) + 0(«-^ 2 ), 

where the "O" term is uniform in x as n — » +oo. □ 

This allows us to blow up the function E n -\(mzw) when z and zv are close to 
the point 1 and wz, n — > +oo with n = m + O(l). 

Lemma 5.3. Ffx a positive real e. For complex rj e C, we ffren /zawe 

E„_i (m(l + m- 1/2 ^)(l + m- ll2 f)j) 

i '- = erf(-£ - fl) + 0(wT2 +e ), 

as m,n —> +oo ty/iiZe n = m + O(l). Here, ffoe "O" expression on the right-hand side is 
uniform on compact subsets ofC. 

Proof. If we put 

C := m(l + m" 1/2 ^)(l + m~ 1/2 fj) = m + m 1/2 (£, + ff) + E,f), 

then 

n-Q n-m- m 1/2 (£, + rj) - £fi ^ / , » 
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where the "O" term is uniform on compact subsets. If we use that E n -\ has only 
nonnegative Taylor coefficients, we get from Lemma 5.2 that 

|E n _!(0| ^ E M _ a (|C|) ___Jn-W 

Vrr 

with a uniform "O" term. So, for real £, r\, we may deduce from | |5.2) the assertion 
of the lemma with e = 0. For general complex '£,, rj, we note that 



(5.2) ^<^=erf(^Wn 
e ig e lcl V \fn / 



Id = m 1(1 + m- 1/2 <£)(l + mr m f\)\ 

= m + m 1/2 Re[£ + ;]] + Re 5 Re rj + \Qm if + |(Im ?]) 2 + 0(m" 1/8 ) 



uniformly in the domain where max{|£|, l?j|} < m 1 ^. As the right-hand side of | |5.2 
is < | for big n, we see that 



E n -i (m(l + m- 1/2 <5)(l + m- 1/2 ^)) 



3 m+ Vm({+r;)+{^ 



3 -(Im £)(Im 17)+ l(Im £) 2 + l(Im r;) 2 +0(nr 1 / 8 ) 

2 



< 2e^ Im ^~ Im ''l 2 



holds in the domain where max{|£|, |rj|} < m 1 ^, provided m is big enough. We need 
to show that he difference 

E„_! (m(l + m- 1/2 4)(l + m- 1/2 /j)) 

Fm,n(Z,n) ■= " " erf(-5 - tf) 

e m+ Vm(£+'j)+£'| 

is of order 0(m _1 ^ 2+£ ) uniformly as r\ remain confined to some compact subset 
of C. We know that F m<n (4;, rj) = 0(m~ 1/2 ) uniformly when £, 1] e IR with confined 
to max{|£|, |rj|} < m 1 ^ 8 . In view of the calculation we just made, we also have a good 
uniform estimate of F m n (£,, 77) when r\ e C withmax{|£|, |rj|} < m 1 / 8 . By employing 
a standard technique involving the subharmonicity of S, i-> log |F„, itI (4, rj)|, and 
certain classical estimates of harmonic measure, we can show that F OT/n (£,rj) = 
0(m e ~ 1 ^ 2 ) holds uniformly when /] e C belong to a compact subset of C, and in 
addition, 77 £ IR Here, e is a positive number which we can get as small as we like. 
A similar argument with rj in place £, worsens the control to F m ,„(£, 77) = 0(m 2£ - 1/2 ), 
but now the control is uniform when both £,, r\ are both complex and confined to 
some compact subset. The proof is complete with e = 2e. a 



The reproducing kernel for a subspace of the Fock space. We shall identify both 
the right-hand and the left-hand side expressions appearing in Lemma [53] with 
reproducing kernels of certain Hilbert spaces of entire functions. This will be the 
case r = of the proposition below. 
Let us agree to identify 

L 2 (1R_) = {/ 6 L 2 (]R) : f(x) = for x > 0}. 

Proposition 5.4. For r = 0, 1, 2, . . ., the function 

(S, rj) h> -^L f ' H r (t + £ - i])H r (t + 77 - l)e-^dt 

is the reproducing kernel for the Hilbert space B,.[L 2 (]R_)] c A 2 r+1 (C). 
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Proof. Let M(£,,i]) = M,,(£) be the reproducing kernel for B,[L 2 (]R_)]. This kernel 
has M (] e B r [L 2 (R_)] and 

— 1— f /(f)H,.( i -, ] -f ? )e t "^" 2 -^ 2 df = B,.[/](, ] ) 
Vr!(27r) 1 / 4 Jr 

= (B,!/]^,)^^, = (^B^M,,])^ 
for all ijeC and all / e L 2 (R_), which allows us to conclude that 

B^lMrjW) = (r!)- 1 / 2 (27i)- 1 / 4 l ] _ K) , 0] (0H,(i - n - fj) e»Hf4' 2 . 
After applying the operator B, to both sides, we see that (cf . Proposition |4.3) 

1 r° 

M n (£) = — — H r (t - £ - De'^-^Hrit -Tj- fj)e tr ^ 2 -^dt 
r\ v2n J-cx, 

= Hr{t-Z-l)H r {t- n -f{)e-^-^&t 

iTl J-ca 

= ^P= H r (f + f7-C)H r (f + ^-^)e^df. 

r! V27T J-oo 



r\ V2t 



The proof is complete. 



Remark 5.5. The special case r = of the kernel in Proposition 5.4 is e 4 ''erf(-£ - 7]), 
which is what we have on the right-hand side in Lemma |5~3] 



The blow-up of the polynomial space at the boundary point. We turn to the 
polyanalytic analogue of the left-hand side of Lemma [B3] 



Definition 5.6. We introduce the blow-up space at 1, 

PoCq := je- ml ' 2 V(l + m~ 1/2 Q : p e Pol m , M ), 
which we equip with the norm 

||| ^ e" ml/2 V(l +^- 1/2 «)|| 6Pol a> ? := m^ 2 eHp\\ L 2 {C ^ r 
For < r < q - 1, we denote by SPoL/ . the subspace 

6Pol iUi := (e" ml ' 2 *P(l + ™" 1/2 £) : p e <5Pol, w+1 }, 
equipped with the same norm. 

An elementary change of variables argument allows us to identify the norm 
on Pol^ M with that of A 2 q (C): 

(5.3) ||| » e-'" 1/2 V(l + ^" 1/2 ^)||p ol a> = \\t » e^Vl + m^)^ (cr 

As a consequence, we may regard PoL^„ . and fiPol^ n r+1 as norm closed subspaces 
of A? (C). We may read off from the definition of the norm in <5Pol a> , , that the 

l,q v ' J m,n,r+x 

kernel on the left-h and side in Lemma 5.3 is the reproducing kernel for the space 

/ -i \ i I 1 1 

Pol .. So, Lemma 5.3 can be understood as saying that 

(5.4) K Pol u> (L n) = K MW] (L q ) + Q(m-i + % 
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where Kp i mA , and ^b [l 2 (r_)] denote the reproducing kernels of the spaces in the 
subscripts, and the bound is locally uniform on compact subsets. We want to 
generalize \5A\ beyond r = 0. To this end, we make use of the operators T r . 



Proposition 5.7. For r = 1, 2, 3, . . ., we have that 

r" 1/2 Ti : SPoli 1 ^ -> <5Pol a> x1 
is an isometric isomorphism. 

Proof. From the isometry properties of T r -\ and T, together with the semi-group 
property r~ ll2 T\ o T r _i = T,-, we get that 

r^Ti: (C) -> ^ (C) 



is an isometric isomorphism. In view of 1 |5.3) , the isometry part of the assertion 
follows. It remains to show that the operator is onto. This is an algebraic exercise 
which we leave to the reader. □ 

By iterating Proposition |5.7{ we obtain the following. 
Corollary 5.8. For r = 1, 2, 3 . . ., we have that 



T r : Pol w , -> 6PoP 

' m,n,\ m,n,r+l 



is an isometric isomorphism 



The blow-up of the polynomial reproducing kernel at a boundary point. From 
Corollary |K8] above, we get that 

(5.5) K 6FolW ($,n) = [T r U%]JK Poim (£,;/)), 

where the subscripts z and w are used to indicate that the operator is acting with 
respect to that variable, and the bar means complex conjugation of the operator. 
To be more precise, 

TV[/](z) := -i=e' z ' 2 ^j/(z)e-' z ' 2 }. 

We would like to plug in the approximation \5A) into | |5.5) . The operator T r is a 
sum of certain polynomials in z of degree < r times the differential operator d z 
to powers < r. The Cauchy integral formula allows us to control the size of the 
derivatives on a compact subset in terms of the size of the functions on a slightly 
bigger compact subset. This means that the approximation < |5.4) carries over, and 
we find that 

(5.6) X 6Pol<1) (L Tj) = 7?)) + 0(m-i + % 

m,n,r+\ * ' 

with uniform control on compact subsets. Next, as 

Poia> m 6Po i<i> 

m,n,q \±s m,n,r+V 
r=0 
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we get that 

9-1 

(5.7) K pol<1 , (^) = YX Pol a> (£,r?) 

1-0 

9-1 
r=0 

again with uniform control on compact subsets. Next, it should be rather clear 
that 

[Tr]{[T r ] r) (K Bo [ i 2 (K _)](^, 7])) 
is the reproducing kernel for the space T r Bo[L 2 (]R_)] = B,.[L 2 (]R_)], which was 



identified in terms of Hermite polynomials back in Proposition 5.4 We write this 
down as a proposition. 

Proposition 5.9. Fix a positive real number e. Then the reproducing kernel for Pol^„„ 
has the following form: 

K Pol u> & 1 \) = Y J —F= H r (t + ^-r])H r (t + i]-Qe-i t2 dt + 0(m-i + % 

asm,n — > +oo while n = m + O(l), where the control is uniform on compact subsets. 
If we like, we may use the classical Christoffel-Darboux identity 

^ L J*'®*™ = ( 9 -l)!(x-y) 

to rewrite the above sum. Also, we should note that reproducing kernel for the 
blow-up space Pol^„ /Cj is connected with the reproducing kernel K M/ „ iCj for Pol m<M 
via the identity 

(5.9) JC Pol<1 » (£, rj) = m- l e- m - mm ^K m , n , q {l + m^, 1 + m^r,). 



The blow-up of the 1-point intensity near a boundary point. The 1-point intensity 
function is 

K (z z1e" m|z|2 

and the localized version with z = 1 + m~ 112 ^ is 

IW3 : = m- x K m ^{\ +m- 1/2 £,l + m- ll2 t)e nil+nrV2 ^, 
where we throw in a factor of m~ x to compensate for the Jacobian. In view of | |5.9) 



together with Proposition 5.9 we obtain 



f -1 n- 2Re£ 

LW# = E — F= H r (f) 2 e-J f2 dt + 0(m-5 +£ ). 

r! V27T J-oo 

So, essentially, the 1-point intensity function is determined by the density 

t k, £ * Hr (*jVi<*, 
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which corresponds to filling the lowest energy eigenstates of the harmonic oscil- 
lator. By the Wigner semi-circle law, then, we get the approximation 



2 r ^Re{ 

— Vl - T dl, 

71 J-l 



valid for big m,n with n = m + O(l), and big q (but much smaller than m,n). So, 
if we rescale to characteristic distance q^^nT 1 ^ 2 we find an interesting law in the 
limit. 

The blow-up Berezin density at a boundary point. The blow-up Berezin density 
at 1 is given by 



1X^(1 + wr 1/2 £,l)l 2 



From 1 5.9 1, we have that 



K Pol fl> (0,0) = m- 1 e- m X m , M (l,l), 



while Proposition |5.9| gives 

*m2> (0,0) = T — = Hrtfe-Pdt + 0{m-\ + % 

as m, n — > +oo with w = m + O(l). Now, as each Hermite polynomial H r is either 
even or odd, 

-0 1 ^+°° r! y^ 



H,.(f) 2 e-2 f2 df = - H r (t) 2 e-i t2 dt 

CO \J — OO 



which leads to 



and 



K Foim (0,0)= q -+O(m-i +£ ) 



K m , n , q (l,l) = | ™?e m [l + 0(m"2 +£ )]. 



1,9^' x ^ ~~ 2 

A similar calculation gives that 
K m , n , q (l +wr x H,\) = me m+mV2 % ^ (60) 

9-1 



1/2 4V — "== f H,.a + £)H,.(t-£)e-5 f2 dt + 0(m-5 +£ )). 



Putting things together, we obtain the following asymptotics for the blow-up 
Berezin density. 



Theorem 5.10. Fix a positive real number e. Then the blow-up Berezin density at 1 has 
the following form: 

9-1 1 f"" 5 2 
Yj~ H,.(f + £)H,.(f-|)e-^ 2 d£ + 0(m^ +£ ), 



1 

Tirj 



asm,n — > +oo zp/zz/e n = m + O(l), w/iere i/ze control is uniform on compact subsets. 



Remark 5.11. When we make some explicit calculations based on Theorem 5.10 



we see that the Fresnel zone pattern is less pronounced for a boundary point z. 
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